Latest supernovae constraints on f{R) cosmologies 
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A class of modified gravity, known as /(7?)-gravity, has presently been applied to Cosmology as 
a realistic alternative to dark energy. In this paper we use the most recent Type-la Supernova 
(SNe la) data, the so-called Union sample of 307 SNe la, to place bounds on a theory of the form 
/(7i) = R — P/R" within the Palatini approach. Given the complementarity of SNe la data with 
other cosmological observables, a joint analysis with measurements of baryon acoustic oscillation 
peak and estimates of the CMB shift parameter is also performed. We show that, for the allowed 
intervals of n, Qmo, and /3, models based on f{R) = R — 13/ R" gravity in the Palatini approach can 
produce the sequence of radiation-dominated, matter-dominated, and accelerating periods without 
need of dark energy. 

PACS numbers: 



I. INTRODUCTION 

One of the key problems at the interface between fun- 
damental physics and cosmology is to understand the 
physical mechanism behind the late-time acceleration of 
the Universe. In principle, this phenomenon may be 
the result of unknown physical processes involving ei- 
ther modifications of gravitation theory or the existence 
of new fields in high energy physics. Although the latter 
route is most commonly used, which gives rise to the idea 
of a dark energy component (see, e.g., [ij), following the 
former, at least two other attractive approaches to this 
problem can be explored. The first one is related to the 
possible existence of extra dimensions, an idea that links 
cosmic acceleration with the hierarchy problem in high 
energy physics, and gives rise to the so-called brane- world 
cosmology [2]. The second one, known as f{R) gravity, 
examine the possibility of modifying Einstein's general 
relativity (GR) by adding terms proportional to powers 
of the Ricci scalar R to the Einstein-Hilbert Lagrangian 
The cosmological interest in f{R) gravity comes 
from the fact that these theories can exhibit naturally an 
accelerating expansion without introducing dark energy. 
However, the freedom in the choice of different functional 
forms of f{R) gives rise to the problem of how to con- 
strain on theoretical and/or observational grounds, the 
many possible f{R) gravity theories. Much efforts within 
the realm, mainly from a theoretical viewpoint, have been 
developed so far [1| (see also Refs. [1] for recent reviews), 
while only recently observational constraints from several 
cosmological data sets have been explored for testing the 
viability of these theories [1, 0, i, [Hl- 
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An important aspect worth emphasizing concerns the 
two different variational approaches that may be followed 
when one works with /(i?) gravity theories, namely, the 
metric and the Palatini formalisms (see, e.g., Q). In the 
metric formalism the connections are assumed to be the 
Christoffel symbols and variation of the action is taken 
with respect to the metric, whereas in the Palatini varia- 
tional approach the metric and the affine connections are 
treated as independent fields and the variation is taken 
with respect to both. In fact, these approaches are equiv- 
alents only in the context of GR, that is, in the case of 
linear Hilbert action; for a general f{R) term in the ac- 
tion they give different equations of motion. 

In the present paper we will restrict ourselves to 
the Palatini formalism for gravitation and will focus on 
its application to the fiat Friedmann-Robertson- Walker 
(FRW) cosmological model. We will derive constraints 
on the two parameters n and /3 of the f{R) = R — /3/i?" 
gravity theory from the most recent compilations of type 
la Supernovae (SNe la) observations, which includes 
the recent large samples from Supernova Legacy Survey 
(SNLS), the ESSENCE Survey, distant SNe la observed 
with HST, and others, giving a sample of 307 SNe la 
events [l2] . We also combine the SNe la data with infor- 
mation from the baryon acoustic oscillation (BAO) [l^ 
and the CMB shift parameter in order to improve 
the SNe la bounds on the free parameters of the theory. 

II. PALATINI APPROACH 

The action that defines an f{R) gravity is given by 

^ " i / '^^^^afiR) + Sm , (1) 

where = SttG, g is the determinant of the metric ten- 
sor and Sm is the standard action for the matter fields. 
Treating the metric and the connection as completely in- 
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dependent fields, variation of this action gives the field 
equations 



(2) 



where T^^^ is the matter energy-momentum tensor which, 
for a perfect-fluid, is given by T^i, — {p„i + Pm)u^u^ + 
Pm5/ii/i where pm is the energy density, Pm is the fluid 
pressure and is the fluid four-velocity. Here, we adopt 
the notation /' = df/dR, f" = d^f/dR'^ and so on. In ^ 
Rfiu is given in the usual way in terms of the independent 
connection rj^^,, and its derivatives, which is related with 
the Christoffel symbol {f^^} of the metric g^^, by 



r^. = {^.} + 2^ (^^9, + d^d^ ~ g^.g'^d,) f (3) 

and R = g^'^Rf,^. 

We assume a homogeneous and isotropic FRW uni- 
verse whose metric is = diag{— 1,0^,0^,0^), where 
a{t) is the cosmological scale factor. The generalized 
Friedmann equation can be written in terms of red- 
shift parameter z = ao/a — 1 and the density parameter 
V,jno = Kpmo/(3-ffo) i^^^ I^^f- for details) 
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where 
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and pmo is the matter density today. The trace of Eq. 
([2) gives another relation 
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and, as can be easily checked, for the Einstein-Hilbert 
Lagrangean (/ — R) Eq. ^ reduces to the known form 
of Friedmann equation. 

By assuming a functional form of the type f{R) = 
R — (3 /i?" , one may easily show that Eq. © evaluated 
at z = imposes the following relation among n, Qmo 
and /? 
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FIG. 1: Hubble diagram for 307 SNe la from the Union sam- 
ple [13] • The curves correspond to the best-fit pairs of n and 
Qmo arising from statistical analyses involving SNe la (only) 
and SNe la + BAG + CMB shift parameter. For the sake of 
comparison the flat ACDM scenario with ^Imo ~ 0.26 is also 
shown. 



III. OBSERVATIONAL ANALYSES 

Since the very first results showing direct evidence for a 
present cosmic acceleration (using a small number of SNe 
la events) [1^1, the number and quality of SNe la data 
available for cosmological studies have increased consid- 
erably due to several observational programs. The most 
up to date set of SNe la has been compiled by Kowal- 
ski et al. and includes recent large samples from 

SNLS [3 and ESSENCE [13 surveys, older data sets 
and the recently extended data set of distant supernovae 
observed with HST. The total compilation, the so-called 
Union sample, amounts to 414 SNe la events, which was 
reduced to 307 data points after selection cuts. 

In this Section, we will use this SNe la sample to place 
limits on the n—flmo (or, equivalently, P—i^mo) paramet- 
ric space. This analysis, therefore, updates the results of 
Refs. [^, [3|. We also perform a joint analysis involving 
the Union SNe la sample and measurements of the bary- 
onic acoustic oscillations (BAO) from SDSS Ts] and the 
CMB shift parameter as given by the WMAP team [lj| 
to break possible degeneracies in the n — ilmo plane (for 
more details on the statistical analyses discussed below 
we refer the reader to Ref. [isj). 



where Rq, the value of the Ricci scalar today, is deter- 
mined from the algebraic equation resulting from equat- 
ing ^ and ^ ioi z — 0. Hence, specifying the values of 
two of these parameters the third is automatically fixed. 
In other words, in the Palatini approach, the two param- 
eter of f{R) — R — (3/R"' can be thought as the pair 
{n,(3) or {n,Vt.-rao)- 



A. Latest SNe la constraints 

The predicted distance modulus for a supernova at red- 
shift z, given a set of parameters P = (n, ilmo), is 

po(z\P) =m- M = b\ogdL + 2b, (8) 
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FIG. 2: Confidence intervals at 68.3%, 95.4% and 99.73% in the n x Qmo plane arising from: Left) SNe la {Union sample); 
Middle) BAG (SDSS); Right) GMB shift parameter (WMAP5). In all panels, the dot marks the best-fit pair (n, fimo) for 
each analysis. Note that the contours of the allowed parameters from these measurements are roughly complementary in the 
n X il,mo space. 



where m and M are, respectively, the apparent and ab- 
solute magnitudes, and di stands for the luminosity dis- 
tance (in units of megaparsecs) , 

where H{z; P) is given by Eqs. ([1]) - (O. 

We estimate the best fit to the set of parameters P by 
using a statistics, with 

_J^ K(z|P)-M^,,(.)]^ 

XsNe — /^ 2 ' y'^^) 

where /^p(z|P) is given by Eq. ([5]), /^^(z) is the extinction 
corrected distance modulus for a given SNe la at z^, and 
Gi is the uncertainty in the individual distance moduli. 
Since we use in our analysis the Union sample (see (T2j 
for details), N = 307. 

Figure (1) shows the Hubble diagram for the 307 SNe 
la events of the Union sample. The curves stand for the 
best- fit f{R) models obtained from SNe la and SNe la -t- 
BAO -1- CMB analysis. For the sake of comparison, the 
standard ACDM model with Q,mo = 0.26 is also shown. 
Note that all models seem to be able to reproduce fairly 
well the SNe la measurements. In Fig. (2a) we show 
the first results of our statistical analyses. Contour plots 
(68.3%, 95.4% and 99.7% c.l.) in the n x Vt„,o plane 
are shown for the given by Eq. p^ . We clearly see 
that SNe la measurements alone do not tightly constrain 
the values of n and ^mo^ allowing for a large interval of 
values for these parameters, with n ranging from -1 to 
even beyond 1, and consistent with both vacuum 
solutions (f^mo = 0), as well with universes with up to 
90% of its energy density in the form of non-rclativistic 



matter. The best-fit values for this analysis are f^mo — 
0.57 and n = 0.99, with the reduced Xr = Xminl^ — 1-01 
{v is defined as degrees of freedom). 

B. Joint analysis 

1. BAO 

The acoustic oscillations of baryons in the primordial 
plasma leave a signature on the correlation function of 
galaxies as observed by Eisenstein et al. IJ]. This sig- 
nature furnishes a standard rule which can be used to 
constrain the following quantity: 

A = Dv^^^^^, (11) 

where the observed value is Aobs = 0.469 ± 0.017, 
z* = 0.35 is the typical redshift of the SDSS sam- 
ple and Dv is the dilation scale, defined as Dy = 

1/3 

[Dm{z*\^)^z^,/H{z^,\P)] with the comoving distance 
Dm given by Dm(z*|P) = J^' dz' / H{z'\P). In Fig. 
(2b) we show the confidence contours (68.3%, 95.4% and 
99.7% c.l.) in the n — flmo plane arising from this mea- 
surement of A. As expected, since this quantity has been 
measured at a specific redshift (z^ = 0.35), it forms bands 
on this parametric space, instead of ellipsoids as in the 
case of SNe la data. 



2. Shift Parameter 

The shift parameter TZ, which determines the whole 
shift of the CMB angular power spectrum, is given by 
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FIG. 3: Left: Confidence intervals at 68.3%, 95.4% and 99.73% in the n x Qmo plane arising from the combined fit involving 
SNe la Union sample, BAO and CMB shift parameter. The best-fit values for this analysis is n = —0.12 and rimo = 0.26 
{(3 = 3.45). Right: Efi^cctive equation of state as a function of redshift for a f{R) = R — 13/R" theory in the Palatini formalism. 
The parameters values correspond to the best-fit solution of our joint analysis. 



— v^^mo 



Hodz' 
H{z';P) 



(12) 



where the z\s — 1089 is the redshift of the last scattering 
surface, and the current estimated value for this quan- 
tity is 7^obs = 1-70 ± 0.03 [3. Note that, to include the 
CMB shift parameter into the analysis, the equations of 
motion must be integrated up to the matter/radiation 
decoupling, z ~ 1089. Since radiation is no longer neg- 
ligible at this redshift, a radiation component with an 
energy density today ofri-y = 5xl0~^ has been included 
in our analysis. Figure (2c) shows the constraints on the 
n — flmo plane from the current WMAP estimate of TZ. 
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3. Results 

In Fig. (3a) we show the results of our joint SNc la 
-f- BAO -I- CMB analysis. Given the complementarity of 
these measurements in the n — Qmo plane [see Figure (2)], 
we obtain a considerable enhancement of the constraining 
power over n and from this combined fit. Note also 
that the best-fit value for the matter density parameter. 



0.26, is consistent with current estimates of 



the contribution of non-relativistic matter to the total 
energy density in the universe (see, e.g., [l3|)- The joint 
fit also constrains the parameters n, flmo, and /3 to lie in 
the following intervals (at 99.7% c.l.) 

n e [-0.3,0.1], e [0.22,0.32] and /? e [1.3, 5.5], 

which is consistent with the results obtained in Refs. [1, 
0| using the supernova Gold and the SNLS data sets, 
respectively. 



TABLE I: Best-fit values for n and /? (the ACDM model cor- 
responds to n = and /3 = 4.38). 



C. Effective equation of state 

Recently, Amendola et al. [20| showed that f{R) de- 
rived cosmologies in the metric formalism cannot produce 
a standard matter-dominated era followed by an acceler- 
ating expansion. To verify if the same undesirable be- 
havior also happens in the Palatini formalism adopted in 
this paper, we first derive the effective equation of state 
(EoS) 



"eff 



= -1 



2{1 + z)dH 

m liz 



(13) 



as a function of the redshift. 

Figure (3b) shows the effective EoS as a function of 
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1 + z for the best-fit solution of our joint SNe la + BAO 
+ CMB analysis. Note that, for this particular combina- 
tion of parameters, the universe goes through the last 
three phases of cosmological evolution, i.e., radiation- 
dominated {w = 1/3), matter-dominated (w = 0) and 
the late time acceleration phase (in this case with w ~ 
— 1). Therefore, the arguments of Rcf. [20] about the 
Weff in the metric approach seem not to apply to the 
Palatini formalism, at least for the interval of parame- 
ters n, Qmo and /3 given by our statistical analysis. In 
Table I we summarize the main results of this paper com- 
pare them with recent determinations of the parameters 
n and /3 from independent analyses. 

IV. CONCLUSIONS 

/(i?)-gravity based cosmology has presently been 
thought of as a realistic alternative to general relativistic 
dark energy models. In this paper, we have worked in the 
context of a f{R) — R — l3R~" gravity with equations of 
motion derived according to the Palatini approach. We 
have performed consistency checks and tested the obser- 
vational viability of these scenarios by using the latest 



sample of SNe la data, the so-called Union sample of 
307 events. Although the current SNe la measurements 
alone cannot constrain significantly the model parame- 
ters n, flmo and /?, when combined with information from 
BAO and CMB shift parameters, the fit leads to very re- 
strictive constraints on the n — rimo (or, equivalently, 
/3 — ilmo) parametric space. At 99.7% c.l., e.g., we have 
found the intervals n G [-0.25, 0.35] and n,no & [0.2, 0.31] 
(/3 e [2.3,7.1]). We note that, differently from results in 
the metric formalism poj , the universe corresponding to 
the best-fit solution for a combined SNe Ia-|-BAO-f CMB 
minimization (ilmo — 0.26 and n — —0.12) shows all 
three last phases of the cosmological evolution: radiation 
era, matter era and a late time cosmic acceleration. 
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